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We exploit theoretical results on the meson spectrum within the framework of a Bethe-Salpeter 
(BS) formalism adjusted for QCD, in order to extract an "experimental" coupling aT P {Q 2 ) below 
1 GeV by comparison with the data. Our results for af^ p (Q 2 ) exhibit a good agreement with the 
£\j ' infrared safe Analytic Perturbation Theory (APT) coupling from 1 GeV down to 200 MeV. 

As a main result, we claim that the combined BS-APT theoretical scheme provides us with a rather 
satisfactory correlated understanding of very high and low energy phenomena. 
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, I. Introduction. The Rcnormalization Group (RG) improved perturbative QCD yields a consistent picture of high 
£SJ ' energy strong interaction processes from a few GeV up to a few hundred GeV scale [1]. At the same time, in the 
, low energy domain the very existence of the unphysical (so-called "Landau" ) singularities in both the RG-invariant 
coupling a s (Q 2 ) and physical observables contradicts the general principles of the local QFT and spoils the theoretical 
analysis of low energy hadron dynamics. In particular, the ghost-pole issue gives rise to severe complications as far as 
O i 1 the bound states problem is concerned, since the scale Q (i.e., the momentum transfer in the qq interaction) involved 
is typically below 1 GeV. Moreover, results of lattice simulations testify to the absence of spurious IR singularities in 
the QCD coupling [2,3]. A reliable algorithm to get rid of these singularities is provided by the APT approach [4], 
based on the causality condition which imposes a s (Q 2 ) to satisfy a dispersion relation with only the physical cut 
1 -oo<Q 2 <0. 

\ This prescription has been exploited in the framework of a second order Bethe-Salpeter (BS) like formalism [5] for 
' the calculation of the meson spectrum in the light and heavy quark sectors. The model is derived from the QCD 
ON ! La grangian taking advantage of a Feynman-Schwingcr representation for the solution of the iterated Dirac equation 
in an external field. Confinement is encoded through an ansatz on the Wilson loop correlator; indeed the quantity 
ilnW^ is written as the sum of a one-gluon exchange (OGE) and an area term 
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>/V i\nW = (i In W) oge +oS. (1) 
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. By means of a three dimensional reduction, the original BS equation takes the form of the eigenvalue equation for 
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M = M£ + UoGE + Ucont , (2) 



where M 2 is the squared bound state mass, M = w\ + W2 = yj m 2 + k 2 + yj m\ + k 2 , k the cm. momentum of 
quarks, mi and m,2 their constituent masses and U — !7oge + C/conf the resulting potential. The combined BS-APT 
theoretical scheme was clearly supported by the results of previous computations performed in [6,7] by using a 1- 
loop APT coupling a^){Q 2 ), with an effective scale constant A^ 1 '^, ~ 200 MeV (see Eq. (4) below). A substantial 
agreement of the spin averaged c.o.g. masses with the data is achieved throughout the whole spectrum and the 
splittings l 3 !?!-! 1 ^ well reproduced in all sectors involving light, strange and charm quarks. Among other attempts 
to study meson properties, by taking relativistic effects into account, we remind e.g. [8] and Refs. therein. They 
differ in the ansatz by which confinement is introduced and in the method used in the determination of bound states: 
quasipotential [8-11], Green function [9] and first order BS formalism [10,11]. 

In this note we summarize the main results of an investigation performed from the reversed point of view, i.e. by 
exploiting the results on the meson spectrum within the BS approach in order to extract "experimental" QCD coupling 
af cp (Q 2 ) below 1 GeV, by comparison with meson mass data. The results are twofold. On the one hand, the 3-loop 
APT coupling reasonably fits a° xp (<3 2 ) from 1 GeV down to 200 MeV, quantitatively confirming the relevance of the 
APT approach to IR phenomena. On the other hand, below this scale, the experimental points give a slight hint 
about the vanishing of a s (Q 2 ) , or the existence of a finite limit lower than l//?o , as Q — > 0. This could correlate with 
some results from lattice simulations [3] and can be theoretically discussed in the framework of a recent "massive" 
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modification [13] of APT. Note that a non vanishing freezing value as suggested in [14] would be still consistent with 
our results for Q < 200 MeV, but does not agree with the general trend of our results in the region 200-500 MeV. For 
the detailed set of our results and technicalities we refer to the extensive account [12]. 

II. The causal APT coupling. A number of non-perturbative tricks to handle the ghost-pole problem was reviewed 
in [15]. Here we exploit the APT approach to QCD (see [4,16]), in which the RG-improved power series in a s (Q 2 ) for 
a given "Euclidean" observable is replaced by a non power expansion over the set of functions 



A n {Q Z ) = 



Si 



da; p n (a) = - lm[a s (-a - ie)] n . 

7T 



(3) 



Here the first function (Xe(Q ) = A\{Q ) plays the role of the APT Euclidean coupling, and at 1-loop it reads 



Po 



1 



A 2 



ln(Q 2 /A 2 ) A 2 -Q 2 



(4) 



At the higher loops Eq. (3) with n = 1 can be integrated only numerically (for details see [17]). Nevertheless, for 
practical applications below 1 GeV one can resort to the same Eq. (4) with modified scale constant (see Rcfs. [12,18]). 
It is relevant to the problem in hand to mention a recently devised "massive" modification for the QCD analytic 
charge [13]. The point is that the representation (3) does not hold for every QCD quantity, and the effect of a 
non vanishing mass threshold m in the dispersion relations could play a substantial role. Then the set of the APT 
functions (3) should be replaced by the set of the "massive" ones with an adjustable parameter m 



A„(Q 2 ,m 2 ) 



Q 2 



Q 2 + Am 2 
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4m 2 



a — Am 2 da 
a + Q 2 ~ 



(5) 



The first function still plays the role of the "massive" coupling a(Q 2 ,m 2 ), with universal limit a(0,m 2 ) = 0. 

III. The quark-antiquark spectrum. Similarly to Rcfs. [6,7], we neglect the spin-orbit and tensorial like terms in 
both the perturbativc and the confining part of BS potential, c^oge and f/conf • Among the spin dependent terms 
only the hyperfine splitting one proportional to |<ti • a 2 is retained. Then one has 



(k|r ()<;E |k') = ^^^p(k,k') 
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(6) 



<k|t/ Conf |k') = 7^3/>(k,k') J d 3 re^- r J inst (r,q, ?10 , g20 ). 
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J inst (r,q, gio ,g 2 o) - 
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Here, q = ^ , Q = k - k' , q j0 = 



(wi+w 2 )(w' 1 +w' 2 ) 



T—r — , while mi and TO2 



- W]+ 2 w :i , Wj = ^Jm 2 + k 2 , and p(k,k') = 

denote the constituent quark and antiquark masses. For the complete expression of the potential and technical details 
we refer to [5,7] and [12]. The meson masses have been computed by the equation 



W 2 = (4> a \M 2 \(/) a ) + ((f>a\UoGE\<i>a) + (<Pa\U C oni\<Pa) , 



(9) 



where (f> a is a zero-order wave function for the state a (a being the whole set of quantum numbers), obtained by 
solving the eigenvalue equation for the static limit Hamiltonian Hqm — wi + w 2 — | + ar by the Raylcigh-Ritz 
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method. To this a second order correction in the hyperfine term was added in some cases. The hurdle of spurious 
singularities has been avoided by replacing a s (Q 2 ) in (6) with a^(Q 2 ) as given by (4) with an effective QCD scale 
A^ 1 '!!^ = 193 MeV. This value has been chosen by imposing that a^(Q 2 ) crosses the 3-loop APT coupling a^(Q 2 ) at 
Q = 0.65 GcV, where a^(Q 2 ) is normalized along with the world average [19] a s (M§) = 0.1176(20) , corresponding 
to A^_ 5 = 236 MeV and A^_ 3 = 417 MeV by continuous threshold matching. The relative difference between the 
1-loop effective and 3-loop exact APT curves is no more than 2% in the region 0.4 < Q < 1.0 GcV, to which the bulk 
of effective Q a belongs, and it is enhanced up to 7% only at Q ~ 0.2 GcV. The string tension has been fixed a priori to 
the value a = 0.18 GeV 2 consistently with lattice simulations. The light and heavy quark masses are then determined 
by fitting the ir , </>, J/tp and T masses [19]. It turns out m u = nrid = 196 MeV, m s = 352 MeV, m c = 1.516 GeV 
and mi, — 4.854 GeV. Within this framework an overall agreement with experimental data is achieved throughout the 
spectrum (see [12] for the complete set of results, preliminary results were given in [20]). 

IV. Extracting a c s xp (Q 2 ) from the data. As stated, we focus our attention here on the reversed problem, i.e., the 
determination of a s at the characteristic scales of a selected number of ground and excited states. As a first step, we 
associate with each state a an effective momentum transfer Q a defined by the relation 

(cf>a\U GE\cf>a) = (4>a \a$ (Q 2 )0(q; Q) \4> a ) = o$ (Q 2 a ) (</>„ |0(q; Q) \4> a ) , (10) 

i.e. as the value of Q for which the fixed coupling value (Q 2 ) inserted in (9) reproduces the same mass m a as 
when using the running coupling a% ] {Q 2 ) . The quantity 0(q; Q) in Eq. (10) can be drawn by the second line of 
Eq. (6). Then, the experimental coupling a c s xp (Q 2 ) can be defined by 

(cb a \M 2 \cj) a ) +aT P (Q 2 a )(MO(q:QM a ) + {^ a \U C onM = m c 2 xp , (11) 
or, by combining Eqs. (9), (10) and (11), 

^M-^M+ Tgkzms - (12) 

The sensitivity of the effective Q's, derived as above, has been checked by analyzing their deviations for a 25% shift 
of A^ 1 ;!^ in (4) and the average change in the momentum scale is about 3%. This shows that a° xp (Qa) is rather 

insensible to the specific form of a^(Q 2 ) , and justifies the use of af^ p (Q 2 ) in (11). Obviously the theoretical meson 
masses are sensitive to a variation of the quark masses (particularly in the case of the ir), while af p and the relative 
Q a turn out to be much more stable. For instance, an increase in the light quark mass of 5% amounts to a change of 
about 2% in the value of af p and 0.2% in the relative Q a ■ 

Note that the APT coupling, involved in the calculation of the spectrum, is remarkably stable with respect to both 
the choice of renormalization scheme and the higher loop corrections [21] (a^(Q 2 ) differs from a^(Q 2 ) by non more 
than 0.3 % below 600 MeV). This makes the method essentially RS independent. 

Our results are unavoidably model dependent due to ansatz (1), which consists of the sum of two contributions that 
one knows to be asymptotically correct for small and large quark-antiquark distances. More sophisticated ansatz also 
exist (see e.g. [22], [23] or [24]), albeit difficult to implement within BS formalism. In the context of our model the 
sources of error are the istantaneous approximation implied by the three-dimensional reduction of the BS equation, the 
approximation introduced into the resolution of the eigenvalue equation, the inclusion of only the leading perturbative 
contribution in the BS kernel /, and finally having neglected coupling between different quark antiquark channels. 
The NLO contribution AI originates essentially from three diagrams with two-gluon exchange; two triangular graphs 
with a four- line vertex g 2 <p*<p A ^A^ and two three-line vertices g<j>*d li <j>A tt , and a crossing box with four three-line 
vertices (Fig. 1). A somewhat crude estimate of these contributions finally yields a global error on the potential 
AO/O ~ AI /I, which spans from 20% for the light-light quark system to about 1% for the bb system. As to the last 
type, within this approximation the BS masses are expected to match the experimental ones within the half width of 
the state, i.e., Am a <~ T a /2 . Keeping in mind Eq. (12), the estimated theoretical errors read 

A NLO a.~aW(^)^ > Ara ^ |(^|0(q;Q)|^)| ra - (13) 

The experimental error Am oxp is generally much smaller than T a /2 . When, however, this is not the case one must 
also consider the experimental error A exp a s , obtained from the second of (13) by replacing m a T a with 2ra e xpAm oxp . 
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Quark self-energy effects have been taken into account by a recursive resolution of the Dyson- Schwinger equation. 
This simply amounts in our approximation to replacing the current quark masses with the constituent masses [7]. 
All other sources of errors above mentioned (including model dependence), though difficult to estimate, can be 
globally taken into account by an additional overall error Am on the masses, independent of a, and chosen such that 
xl = jvSF ^a=i( m a - m c ^ p ) 2 /[(A tot m a ) 2 + (Am) 2 ] ~ 1 . Here A to tTOa is the total error resulting from all sources 
explicitly evaluated, (Atot^a) 2 = m 2 AI/I + (T a /2) 2 + (Am cxp ) 2 and the sum restricted to the safer S and P states. 
We find Am ~ 20 MeV. 

V. Conclusions. All results are displayed pictorially in Fig. 2. Values of a| xp at the same Q from triplet and singlet 
states have been combined by means of a weighted average according to their errors. The points a e s xp (Q 2 ) show 
a noticeable evolution from 500 down to 200 MeV, where only the safer S and P states are involved, in remarkable 
agreement with the 3-loop APT coupling a$ (Q 2 ) properly normalized. Precisely we find x 2 = a^tj J2a=i ( a s Xp (Qa)~ 

ct~E (Q 2 )) 2 /[(A tot Q! s ) 2 + (Aa s ) 2 ] <~ 0.8, with A tot Q! s the total error explicitly evaluated and A« s the uncertainty 
obtained from the second of (13) by replacing T a by 2Am. The agreement quantitatively supports the APT approach 
to IR phenomena down to a few hundred MeV. Below 200 MeV, the experimental points exhibit a tendency to deviate 
from the APT curve and to approach zero or at least a finite limit lower the universal APT freezing value. Note that 
these points have been obtained from high orbital excitations (D and F states), both experimentally and theoretically 
much more uncertain, as clearly shown by their large error bars. However, the extracted points a| xp (Q 2 ) could 
correlate with some lattice results [3], and discussed in the framework of the "massive" modification [13] of the APT 
algorithm (see Sec. 1) which takes into account effects of a finite threshold. Finally we stress that a synthesis of results 
for a s (Q 2 ) as defined from bound states in the BS framework with high energy data shows a very good agreement 
with the 3-loop APT coupling. The pcrturbativc 3-loop coupling with IR singular behavior is ruled out by the data, 
whereas the BS-APT theoretical scheme allows a rather satisfactory correlated understanding of very high and low 
energy phenomena. 
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FIG. 2. Extracted values of aT P (Q 2 ) against the 3-loop APT coupling (3) with > =3 = (417 ± 42) MeV (solid), and its perturbative 
counterpart (dot-dashed). The "massive" 1-loop APT coupling (n = 1 in (5)) refers to A^'Z.g = 204 MeV and m cff = (38 ± 10) MeV 
(dashed). Circles, stars and squares refer respectively to qq , ss and qs with q = u,d, diamonds and crosses to cc and bb; asterisks stay 
for qc and qb , whereas plus signs for sc and sb . 



